Quantum phase transitions in trimerized zig-zag spin ladders 
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We analyze the effects of a trimerized modulation in a quantum spin S — \ zig-zag ladder at the 
magnetization plateau M = f/3. Such periodicity is argued to be stemmed from lattice deforma- 
tions by phonons. The interplay between frustration and exchange modulation is well described by 
an effective triple sine-Gordon field theory close to the homogeneous ladder and by block-spin per- 
turbation theory in the weakly coupled trimers regime. The characteristic triple degeneracy of the 
ground state for homogeneous ladders gives place to modulation driven quantum phase transitions, 
leading to a rich phase diagram including up-up-down, quantum plateau and gapless plateau states. 

PACS numbers: 75.10.Jm, 73.43.Nq, 75.30.-m 



I. INTRODUCTION 



Frustrated quasi one dimensional antiferromagnetic 



spin S 



systems have been extensively studied in 



the last years. One of the most discussed and paradig- 
matic models is the 3\ — J-2 zig-zag ladder, which apart 
from the theoretical interest caused by the frustration is 
believed to be the relevant starting point for describing 
magnetic excitations of a number of real quasi one di- 
mensional materials 1 -. Among others, attention is focused 
on CuGeOs where, along with frustration, spontaneous 
spin-Peierls exchange modulation caused by the coupling 
between spins and phonons plays a crucial role^. Due to 
the interest in the above material the interplay between 
exchange modulation and frustration has stimulated con- 
siderable efforts in purely spin systems as well, with stud- 
ies concentrating on the zero magnetization case2^£&£. 

Recently it was shown that phonons in frustrated zig- 
zag ladders can open plateaus not only at zero mag- 
netization but at other rational values M of saturation 
(M=l/3, l/2...)&. Spontaneous exchange modulations in 
such situations could also take place, with a spatial pat- 
tern associated to the (generally broken) ground state 
translational symmetry. 

In this work we will concentrate on the interplay be- 
tween exchange modulation and frustration in a pure zig- 
zag spin ladder at the M=l/3 plateau state. This system, 
in absence of modulation, is known to exhibit a magneti- 
zation plateau when the next nearest neighbour (NNN) 
coupling Ji is large enough with respect to the nearest 



neighbour (NN) coupling J 



9.10 



Moreover, the ground 



state is three-fold degenerate and translation symmetry 
by one lattice spacing is spontaneously broken to an up- 
up-down configuration 9 -. A natural elastic deformation is 
then given by a period three pattern 



o sm( — % 

y 3 
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where u% is a relevant scalar coordinate describing the 
displacement of the ith ion sequentially numbered on a 
one dimensional chain. We will pay attention to the most 
symmetric situation J\ = J2, and a deformation with 
amplitude 5 > and phase $0 = — -|, that causes site 
1 and 3 to move closer to site 2 and so on. We assume 
that magnetic exchange couplings vary linearly with the 
relative displacement to estimate the modulation on both 
NN and NNN couplings. In this way, sites 1,2,3 are 
magnetically coupled by equally enhanced exchanges J = 
J\ + A forming equilateral triangular trimers as shown in 
Fig. [TJ Couplings drawn with dashed lines are weaker 




FIG. 1: 
ladder, 
lated. 



Schematic description of the trimerized zig-zag spin 
Both nearest and next nearest couplings are modu- 



than Ji, which for simplicity we represent as J' = J\ — A. 

Notice that this deformation selects triangles with 
stronger couplings as basic frustrated units, with weaker 
couplings between them; in a limiting situation, that we 
will also consider, the system is formed by weakly cou- 
pled triangular trimers where strong frustration enhances 
plateau formation and directly relates to high degeneracy 
of the ground state. On the other hand, a deformation 
with 8 < (stronger couplings along dashed lines in Fig. 
[J) modifies the system towards a single homogeneous spin 
chain with modulated weaker couplings between non NN, 
up to fourth neighbours; the corresponding limiting situ- 
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ation (J/ J' = 0) is just a usual spin chain running along 
the dashed line, where no plateau would be observed in 
the magnetization curve. 

Guided by this quick analysis, in the present paper we 
will investigate the M=l/3 ground state of the trimer- 
ized antiferromagnetic spin S = h zig-zag ladder. Specif- 
ically, we start with a homogeneous Hcisenbcrg antifer- 
romagnetic zig-zag ladder with exchanges J\ = J2 at 
1 /3 magnetization. This situation is well above the crit- 
ical coupling J2 = 0.487Ji— , so that the system exhibits 
the magnetization plateau with a three-fold degenerate 
ground state. We then consider a lattice deformation 
of period three, giving rise to a modulation of the same 
order on NN and NNN exchange couplings. The Hamil- 
tonian can be written as 



AL 



H 



(Ji Si ■ iSj+i + Ji Si ■ S'i+2), 



(2) 
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FIG. 2: Schematic description of the ground state phase dia- 
gram of the spin system shown in Fig. [1] M 3 is a magnetiza- 
tion order parameter to be referred to in Eq. (|22[) . 



where Si denotes a spin | operator at site i. The modu- 
lation is given by NN antiferromagnetic couplings Ji > 
forming a sequence of period three with J\ = J 2 = J, 
J3 = J', and NNN antiferromagnetic couplings J; > 
also forming a sequence of period three with J\ = J, 
J2 = J3 = J'- 

As we will see, the interplay between frustration and 
modulation gives rise to a very rich ground state phase 
diagram involving a number of quantum phase transi- 
tions of different order, including up-up-down, quantum 
plateau and gapless chiral plateau states. In summary, 
our analysis will show that the system adopts a unique 
up-up-down ground state for a = J'/ J > 1 separated by 
a first order transition at 0:3 = 1 from a Z2 degenerate 
ground state phase extending to a slightly less than 1. At 
some finite distance «2 < 1 we find a second order tran- 
sition of Ising class to a unique quantum plateau state. 
These quantum phase transitions should be contrasted 
with the stability of the three-fold degenerate ground 
state against translation invariant modifications of the 
zig-zag couplings J\ , Jj^. Towards the strong frustration 
regime a< 1 we find another transition at < a>i < oli 
to a gapless chiral phase with uniform magnetization 
ground state. Thus we provide a phase diagram with 
J' starting from zero (decoupled triangular trimers) up 
to J' > J, schematically shown in Fig. [2] The anal- 
ysis relies on bosonization, renormalization group and 
block-spin perturbations, supported by numerical diago- 
nalization on finite size systems. A detailed knowledge of 
these transitions will ultimately help understanding the 
possibility of spin-Peierls-like phenomenon in magnetized 
materials. 

The paper is organized as follows. In section II we 
study the regime of weak modulation J' ~ J by bosoniza- 
tion, obtaining a triple sine-Gordon effective theory. In 
section III we study the strong frustration regime, where 
block-spin perturbation theory is applicable. Effective 
Hamiltonians at first and second order are analyzed. In 
section IV we study numerically the ground state and 
low lying excitations in finite systems by exact diagonal- 



ization; this supports the semi-quantitative bosonization 
results and provides a bridge between weak modulation 
and strong frustration regimes. In section V we present 
our conclusions and prospects for future work. 



II. WEAKLY MODULATED SYSTEMS 

We first analyze the stability of the zig-zag ladder 
plateau ground state at M = 1/3 against small mod- 
ulated perturbations, as defined in Eq. © for J' ~ J. 
The microscopic Hamiltonian for the system can be con- 
veniently written as 



H = 



E 



[- 



2tt s 



(J - e)Si ■ S t+2 - hSi 



-e^cos{i — )Si ■ S l+ i 

i 

+£^cos((i- l)y)£ 



(3) 



where e = 2(J — J')/3. The first line describes a homo- 
geneous zig-zag ladder with Ji = J — e/2 and J 2 = J — e 
in an external magnetic field h, while the rest is a mod- 
ulated perturbation of period three. 

Once the magnetic field is tuned to the plateau range, 
the homogeneous part is well describe d 10 ' 12 by bosoniza- 
tion with an effective massive sine-Gordon Hamiltonian 



d.r 



K (d x <f>f 



gcos[3v47T0 



. (4) 

where is a real bosonic field with spin wave velocity v 
and compactification radius R = \j\f^^. and 4> is its 
dual field defined by d x (j> = dt<j>. The Luttinger param- 
eter K takes into account renormalization due to spin 
interactions. The presence of the third harmonic of the 
bosonic field arises from a triple Umklapp term, only al- 
lowed at the Fermi momentum kp = 7r/3 fixed by the 
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magnetization M=l/3. Moreover, K has to be suffi- 
ciently small. K < 2/9, rendering this term a relevant 
conformal perturbation, and g has to be positive, in or- 
der to fit numerical results. 

Indeed, in a semiclassical analysis, the relevant third 
harmonic — g cos[3^/iTT(f>] is considered as a potential; it 
has three non equivalent minima for the compactificd bo- 
son field <f> (see Fig. [3]) signaling a three-fold degenerate 
ground state. The mapping to spin variables, following 




FIG. 3: Semiclassical potential for the homogeneous system 
J — J' (in arbitrary units). There are three minima in the 

2 ' 2 > 

non-equivalent up-up-down states. 



compactified range (— ^y-,-*y-), corresponding to the three 



the usual rule for magnetized systems^ 3 -, shows that each 
of the minima for <p corresponds to the so called "classi- 
cal plateau states"—, namely up-up-down states related 
by lattice translations i — > i + 1, i + 2. This is just what 
numerical results^ show, a three-fold degenerate ground 
state with spontaneous Z3 symmetry breaking to states 
with up-up-down local magnetization structure. The de- 
nomination of "classical plateau states" indicates that 
these are essentially the states obtained in the Ising limit 
of the present spin system. Interestingly, this description 
is stable for a wide range of homogeneous variations of 
J 1 and J 2 couplings. 

The second and third lines in Eq. (|3|) represent a mod- 
ulated perturbation to the homogeneous zig-zag ladder 
described above. The key point in our presentation 
is that, after bosonization, these perturbations provide 
first and second harmonics in the Hamiltonian that are 
commensurate with the short scale oscillations given by 
kp = 7r/3. Specifically, these terms yield an effective 
contribution which can be recasted as 



Hmod ~ e J dx C cos[V^7r</>] — cos[2-\/47r(/)] 



(5) 



where C is a coefficient of order 3 (C = 1 + =^-). These 
harmonics have smaller scaling dimension than the third 
one, thus providing relevant perturbations that will com- 
pete with it. The presence of the first harmonic substan- 
tially modifies the theory It is now the leading pertur- 
bation, and the effective field theory Hq + H mo d is the 



so called triple sine-Gordon modet^. Extensive analysis 
of competition between harmonics has been presented 
in [l^fl6llT7| . mainly focused on the double sine-Gordon 
model. The three-frequency case has been recently dis- 
cussed in detail in [l8[ . 

Regarding spin systems, a model similar to ours has 
been recently analyzed by Hida and Affleck in [l2j , where 
a modulated perturbation on NN couplings of a zig-zag 
ladder was proposed as a way to diminish frustration and 
enforce singlet formation on the M=l/3 plateau ground 
state. The effective theory obtained by these authors 
contains first and third harmonics of the bosonic field. In 
the present case, the second harmonics makes necessary 
an extension of the above mentioned results. 

We will first perform the semiclassical analysis of the 
ground state, using the bare coefficients in Eq. (J5|). 
Thereafter, the Rcnormalization Group (RG) flow of the 
relevant perturbations will be discussed. 

For J 1 > J the basic harmonic dominates the mini- 
mum of the potential (see Fig. [H upper panel), so that 
the perturbation selects only one configuration = 
corresponding to one of the up-up-down states out of the 
three-fold degenerate ground state (the one with spins 
down on sites numbered by 3i + 2 in Fig. [1} . 

A richer scenario occurs for J' < J when the first and 
second harmonics are in conflict with the triple harmonic 
term (see Fig. |4j lower panel) . In this case triple degen- 
eracy is transformed first into double by an infinitesimal 
perturbation. From the three degenerate states at J' = J 
it is precisely the state selected for J' > J the one that 
now becomes excited, while the other two potential min- 
ima are selected and shifted from their commensurate 
up-up-down positions. This phase is characterized by Z2 
reflection symmetry, spontaneously broken by each of the 
potential minima. When J — J' reaches some finite pos- 
itive value these two minima collapse into a single one at 
4> = = — 0T/2, lifting the degeneracy completely. 

Moreover, the position of the minimum is shifted to a 
commensurate field configuration which was a maximum 
for the homogeneous ladder. In terms of spins, the system 
is pinned in a very different configuration, the so called 
"quantum plateau state"—. This configuration is charac- 
terized by spin singlets alternating with spin up sites (in 
our case, with spin up at sites numbered by 3i + 2 on Fig. 
[1]). Given the symmetries and degeneracy of the ground 
states on both sides of the critical point, one may conjec- 
ture that the transition belongs to the second order Ising 
universality classic. 

The preceding discussion describes the competition be- 
tween harmonics, depending on the signs and values of 
their bare coefficients. The results are particularly sen- 
sitive to the relative coefficients of first and second har- 
monics. In order to describe the long distance effective 
theory, the RG flow of the coefficients must be analyzed. 
To this aim, let us write the triple sine-Gordon Hamilto- 
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and g3, the perturbative RG equations read 



J'/J > 1 
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FIG. 4: Modification of the semiclassical potential by relevant 
harmonics (in arbitrary units) . For J' > J (upper panel) , the 
single central minimum pins the system in a particular up-up- 
down state. For J' < J (lower panel), the minima structure 
takes the ground state first to double degeneracy and then, at 
some finite critical value of J'/J, to a single minimum shifted 
\pK 1% with respect to non-perturbed position. This shift se- 
lects a quantum plateau state out of the triple degenerate 
up-up-down states. 



nian in a compact form 



H = Hll + <7i cos[v47T0] — g 2 cos[2v Att4>] — 

—g 3 cos[3"\/47r</>] (6) 



dlK 

dm 

dl 
dgi 
dl 

dm 
dl 
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(2 - K )gi + — g x g 2 - ^5233 
4 

(2 - AK)g 2 - - gi g 3 , 



(7) 



where I is the length scale. The first two equations arc the 
well-known sine-Gordon RG equations^; as expected, 
<?3 flows towards strong coupling and K decreases. The 
other two equations can then be solved using the solu- 
tions for K,g%. We performed a numerical analysis of 
Eqs. in the region of interest. The typical form of 
the RG flow is shown in Fig. [5l The salient features are 
that <7i remains positive and grows faster than g$, al- 
lowing for the modification of the semiclassical potential, 
and that in the perturbative regime 172 remains positive 
and smaller than g\ . This analysis proves that the phase 
transitions depicted above are correct in the low energy 
(long distance) renormalized effective theory. 
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where Ull is the free boson Hamiltonian (Luttingcr liq- 
uid) with Luttingcr parameter K . All the present cosine 
terms are relevant perturbations when K < |, and the 
signs have been chosen so as to represent the bare sit- 
uation with positive coupling constants 31, 32, 33- The 
region of parameters we are interested in corresponds to 
#3 of order of unity, describing the homogeneous ladder 
Z 3 symmetric ground state, and small g±, g 2 oc e describ- 
ing the modulation effects. Up to second order in g\, g 2 



FIG. 5: Representative RG flow of triple sine-Gordon cou- 
plings in terms of the length scale. The region of interest is 
given by initial conditions K < | , g% 3> <?2 , fli and gi ~ 3g2 , 
corresponding to the bare Hamiltonian in Eqs. (T3] [5J). (a) Cou- 
plings K and gs describing the homogeneous zig-zag ladder, 
(b) First and second harmonics couplings gi and g^ associated 
to modulated perturbations. 

We have then found that the bosonization analysis pre- 
dicts a first order phase transition at J' = </, as the 
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configurations selected for ground states by the classical 
potential jump between different minima of the poten- 
tial, and a second order phase transition of Ising class 
at some finite value of «2 = J' /J, where two degenerate 
minima merge continuously into a single one. 

A natural question is now if the quantum plateau phase 
remains stable for J'/ J < oci- It is easy to see that there 
should be a further phase transition, as for J'/ J <C 1 
the system is described by almost decoupled triangular 
trimers. We investigate this phase in the following sec- 
tion. 



III. STRONG FRUSTRATION 

In order to analyze the ground state of the system be- 
yond the weak modulation regime, we choose in this sec- 
tion a different starting point, namely the strong frus- 
tration regime. We consider the system as consisting of 
frustrated triangles of spin S = \ with antiferromagnetic 
exchange coupling J, weakly coupled among them with 
exchange coupling J' <C J. In this regime, the system 



1 J' 2 3 




2 J 3 1 



FIG. 6: Schematic description of the spin system in the strong 
frustration regime, with J' < J indicated by bold line 
trimers. 

can be thought of as a quasi one-dimensional weakly cou- 
pled chain of triangles. The topology of weak couplings 
has been chosen so as to recover the zig-zag chain in the 
limit J' I J — * 1, but similar systems with more symmetric 
inter-triangle couplings have been considered in relation 
with spin tubes2£. Notice that the M=l/3 plateau is very 
robust for J 1 <C J, as becomes apparent in the limit case 
J' = 0; we can then safely consider a plateau regime with 
M=l/3 even for small magnetic fields. 

For convenience, in this section we enumerate the tri- 
angular trimers with an index n and rename spins as S% 
(a = 1,2,3 inside each trimer, see Fig. O to write the 
Hamiltonian in Eq. ([2]) as 

H = ^ J (^S n ■ S 2 + ■ + S„ ■ §n\ (8) 

n 
n 

We perform a systematic block-spin perturbative 
analysi o 21 i 22 around the highly degenerate exact ground 
state of the system of decoupled trimers ( J' = 0) at zero 
magnetic field. For low energy, the spin operators at each 



vertex can be factorized in terms of the triangle total spin 
operator and a pseudo-spin chiral operator— as 

Sn = \Su,t{\ ~ T%), (9) 

where 

S n ,T = S n + S n + S n (10) 

is the total spin operator of the n-th triangle projected 
onto the S = \ low energy sector. Operators act on 
a chiral sector as 

T 1 = T++T- 
n n 1 n > 

T n = U 2 T+ + UJT~ , 

Tl = UJT++^T-, (11) 

where r„ are generators of the 5* = \ pseudo-spin sector, 
with t* = -^S n ■ S n A S n , and u> = e l2jr/3 a primitive 
cubic root of unity. 

In order to study the M=l/3 plateau ground state, 
we consider the system under the action of an external 
magnetic field h* high enough so as to force total magne- 
tization 1/3, but still low enough to discard excited spin 
states S = | at any triangle. This regime is clearly avail- 
able for small J' / J. The trimer spin degree of freedom is 
then saturated to T = \ , leaving a pseudo-spin chain 
Hamiltonian describing non magnetic excitations. Notice 
that the magnetic field couples only to the total spin and 
plays no role in the pseudo-spin sector. 

At first order in block-spin perturbations we get a sim- 
ple XY-like nearest trimers Hamiltonian, 

H Sf = 1T J 'E I 6 " 1 T n Vfl + <^T+ +1 T-), (12) 

n 

where we have dropped constant terms, including the 
magnetic field. One can perform a gauge transformation 

< - r+, 

< - <, (13) 

that preserves SU(2) commutation relations, to eliminate 
phases. Then the usual Jordan- Wigner transformation 
maps non magnetic chiral excitations onto a gapless the- 
ory of free spinless fermions. The ground state is thus 
described by a half-filled fcrmion band. 

From this effective model, one can also analyze the 
non magnetic excitations above the plateau ground state, 
which form a gapless continuum with non degenerate 
ground state. Using Eqs. (|10ll 1I13|) one can evaluate 
the ground state expectation value and correlations of 
spin operators. The spin density in the direction parallel 
to the applied field is uniform and simply equals average 
magnetization, 

= 1/6, (14) 
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whereas in-plane averages vanish. Equal time spin-spin 
correlation function in the direction parallel to the ap- 
plied field behaves as: 



fectivc Hamiltonian written as 



/ qz a qz b\ 



const 
^/n — m 



2k 



cos((n - m)- + —(a - &)), 
(15) 

with a, b = 1,2,3 and thus show long range order (each 
triangle being in a spin +1/2 state), plus algebraically 
decaying oscillations. The in-plane- ATY correlation func- 
tions decay exponentially. This picture is a consequence 
of the fact that non magnetic excitations are gapless, 
while magnetic excitations are gapped. 

From the above discussion one can state that for J' <C 
J the spin system at M = 1/3 presents a gapless phase 
corresponding to non magnetic chiral degrees of freedom 
described by a Luttinger liquid with K = 1. 

In order to look for a quantum phase transition to- 
wards the quantum plateau state, it is necessary to con- 
struct the block-spin effective Hamiltonian at second or- 
der. This derivation requires a much longer calculation, 
and provides next nearest trimers interactions. Skipping 
details, once total spin is saturated to \ at each triangle 
we get for the second order correction 



H eff/ J 



< 27 X^ 7 ™ 



(16) 
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where again we have dropped constant contributions. 
Terms in the second line are similar to those appear- 
ing at first order in Eq. (fT2)) . so we propose a different 
gauge transformation in order to eliminate the phase in 
XY terms 



T-+ _j p iX ( a ) n t+ 

'n ° ' n ' 



(17) 



Heff/J = 

-P(«) Y tj( T « T "+l + T n+l^n) + A 2 (a)« +1 ] 



-i2(X 



{a) ~* ] ^Tn+2+h.c:_ 



81 



a 2 ^ [e i3 A (a) („ +1 ) r+T + +iT++2 + ^ (18) 



where A z (a) 



3p(a) 



and p(a) = a/8lV972 + 25a 2 



In Eq. (fl8|) the second and third lines represent an ex- 
tended Heisenberg XX Z model with NN and NNN inter- 
actions for pseudo-spin operators. NN interactions have 
acquired the anisotropy parameter A z (a), while NNN in- 
teractions are XY-like, with a shifted band (because of 
the phase in the couplings); inspection of the coefficients 
shows that < A z (a) < 1 and NNN couplings are small 
compared with NN ones. All the other terms include 
phases that depend on the position, and will generally 
cancel out in the continuum limit, so that the continuous 
U(l) symmetry broken by some of them is recovered. 
The low energy behaviour of this pseudo-spin i model 
is then appropriately described by bosonization, where 
it is most clearly seen that oscillatory terms are incom- 
mensurate and can be neglected. We arc essentially left 
with a one-dimensional spin S = | anisotropic Heisen- 
berg model. 

Applying the standard bosonization procedure we de- 
rive the low energy Hamiltonian 



eff 



\ J dx[±(d x <p) 2 +K(d xl p) 
-T J dx cos[2\/4mip\, 



(19) 



where the Fermi velocity, the Luttinger parameter and 
the coupling T depend on J' /J. For small J'/J the har- 
monic perturbation is strongly irrelevant, so that the ef- 
fective theory still describes a gapless phase. As J'/J 
increases, the conformal dimension of the harmonic per- 
turbation decreases and one expects from Eq. (|19jl a 
second order Brezinskii-Kosterlitz-Thouless (BKT) like 
phase transition to a massive phase. A detailed compu- 
tation allows to estimate a critical value J'/J ~ 0.5, a 
value that is possible beyond the validity of block-spin 
perturbation theory. 



IV. NUMERICAL ANALYSIS 



where X(a) = arctan[ 5 l ^ ] , rendering the total ef- 



In order to support the analytical results in the pre- 
ceding sections, and to explore the intermediate regime 
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of J' / J not covered by perturbations around J' = and 
J' = J, wc performed exact diagonalization analysis on 
finite systems with 12, 18 and 24 spins using periodic 
boundary conditions. However, in attempting to avoid 
misleading results in the approximants and extrapola- 
tions referred to below, we discarded the 12 spin data. 

We have first confirmed the presence of the M = 1/3 
magnetization plateau, all over the range from J' = to 
J' = 2 J. The magnetic phase diagram showing the mag- 
netic field h necessary for level crossing between available 
magnetizations in several finite size systems is presented 
in Fig. [7J A noticeable finite increment in h separates 
magnetic excitations from the M = 1/3 ground state, in- 
dicating the magnetization plateau. This is interpreted 
to remain in the infinite size scaling limit. 




FIG. 7: Magnetic phase diagram obtained by exact diagonal- 
ization of finite system with 12 (doted line), 18 (dashed line) 
and 24 spins (solid line). Inset: a typical magnetization curve. 

We have then computed the ground state and first 
three excitation energies in the subspace of magnetization 
M = 1/3 in a wide range of couplings < J'/ J < 1.5. 
In Fig. [8] we plot the gaps to excited energies in terms of 
J' I J. The triple degeneracy of the zig-zag ladder ground 
state is, within finite size effects, qualitatively observed 
at point J' I J = 1 in agreement with Ref. 0. More- 
over, assuming that the ground state will become degen- 
erate in the thermodynamic limit, our numerical data are 
compatible with the picture that the triple degeneracy is 
lifted to a unique ground state for J'/ J > 1. It also seems 
to be partially lifted to double degeneracy for J' / J < 1 
where one level rapidly separates but another gets closer 
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FIG. 8: Non magnetic excitations X n for 24 spins with pe- 
riodic boundary conditions above the M = 1/3 ground state 
energy Ao, for n = 1 (dotted line), 2 (solid line), and 3 (dashed 
line). The inset shows the corresponding results for 18 spins. 
Numerical data in all subsequent figures were also obtained 
fixing M = 1/3. 



to the ground state, and only then raises lifting the re- 
maining degeneracy. 

One can estimate the locations of the Ising and BKT 
transition critical points mentioned in sections II and 
III by considering the Callan-Symanzik /^-function de- 
veloped in the context of phenomenological renormaliza- 
tion group by Roomany and Wyld [24 |. This technique 
can handle situations in which the phase transition is 
not necessarily characterized by a power law decay of 
the spectrum gap, i.e. an ordinary second order transi- 
tion, but also by a singular BKT form. In the former 
case the /3-function exhibits a simple zero whereas in the 
latter situation it vanishes with an algebraic singularity. 
This function can be estimated from finite lattice data 
by the Roomany- Wyld (RW) approximant, which in our 
notation reads 



0Rw{a) = a In 
'N + 6' 



N + 6 A 



N+6 



(a) 



N 



A N (a) 



/ 



(20) 



In 



N 



1 



ad a ln[(N + 6)A N+6 (a)NA N (a)\ 



where Ajv is the spectrum gap per spin. Notice that 
whenever the phenomenological renormalization condi- 
tion (N + 6)AN +6 (a) = NAjy(a) is satisfied, the [5- 
function shows a zero. Furthermore, it behaves as 
f3{a>) ~ (a — a c )/u, from where the slope at a c is re- 
lated to the exponent of the second order transition. In- 
stead, in the vicinity of a singular transition of the form 
A oc cxp— (a — a c ) _<T , one has (3 ~ (a — a c ) 1+a which 



8 



ultimately determines both a c and a. In Fig. [9] we plot 
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FIG. 9: The Roomany - Wyld approximant (N = 24, 18) of 
the Callan-Symanzik j3- function (circles). The left and right 
curves are respectively consistent with a BKT phase transi- 
tion at a\ ~ 0.35 and with an Ising transition at 012 ~ 0.88. 
Dashed curves fit the numerical data with the parameters re- 
ferred to in the text. 

the RW approximant computed from 18 and 24 spins 
data. The leftmost curve characterizes a singular tran- 
sition with ai ~ 0.35 and a ~ 0.7 (though not exhibit- 
ing a strict zero, possibly an artifact of our small lattice 
sizes) compatible with the BKT type, while the right- 
most one typifies a conventional transition at ~ 0.88 
with v ~ 0.67 consistent with the Ising class. 

The BKT transition can also be estimated from level 



crossing spectroscopy 



25.20.27 



of the low-lying states with 



different symmetries. In Fig. [101 WG show the size-scaled 
excitation energies NAE(N) as functions of J'/J for 
finite-size clusters with N = 18 and N — 24 spins. The 
intersection between the first and second excitations is 
interpreted^ as the chiral fluid-quantum plateau tran- 
sition critical point a±(N) and occurs at a = 0.29 and 
a = 0.36 respectively. The finite size scaling of ai(N) is 
expected to follow 



J25 



cti(N) = ai(oo) + const x N 



(21) 



suggesting a crude extrapolation to ai(oo) w 0.47. 

These numerical estimates are thus consistent with the 
existence of a BKT transition at ct\ in the range 0.3 — 0.5 
and an Ising transition at CV2 around 0.9. However, a 
level crossing, not predicted by the analytical treatment, 
is seen in Fig.[5]at J'/J m 0.18 in the 24 spins system but 
not in the 18 spins one. Its presence should be checked 
in larger systems, not currently available to us, as it may 
well be the consequence of highly oscillating terms in Eq. 
(fT8"]) which become particularly dominant for small lat- 
tice sizes. If confirmed, it would indicate that the BKT 
transition separating the gaplcss phase described in the 
strong frustration regime from the quantum plateau de- 
scribed in the weak modulation regime could be replaced 
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FIG. 10: Detail of scaled excited levels crossing used in the 
level spectroscopy analysis. 



by a first order one in the thermodynamic limit. We can 
not make a definite distinction from our current numeri- 
cal data. 

In order to characterize the phases separated by the 
mentioned transition points, we also computed the lo- 
cal magnetization profile for the ground state. We found 
three different periodic phases according to the generic 
J'/J values explored. As the profile is periodic, we re- 
port the configuration of a generic trimer, labelling sites 
a — 1,2,3 in accordance with Fig. [51 A plot with local 
magnetization of the ground state for J'/J up to 1.4 is 
shown in Fig. [TT] (upper panel). Notice that a = 2, 3 sites 
show the same magnetization. 

In the strong frustration limit we observe for the 
24 spins system an almost uniform magnetized ground 
state for < J'/J < 0.18 (for instance S zl - 2 ' 3 = 
0.018, 0.016, 0.016 at J'/J = 0.1). In contrast, at 
intermediate regimes we find a quantum plateau like 
magnetization for 0.18 < J'/J < 1 (c.f. S zl < 2 ' 3 = 
0.470, 0.015, 0.015 at J'/J = 0.5), as well as an up- 
up-down magnetization for modulated regions J'/J > 1 
(c.f. S zl ' 2 ' 3 = -0.075, 0.288, 0.288 at J'/J = 1.2). 
We find clear signals of a level crossing at J'/J = 1: 
the ground state magnetization profiles correspond to 
a quantum plateau (S z f=a 0,1/2,0 at each trimer) for 
0.18 < J'/J < 1, and up-up-down states (sign(S' 2 ) = 
-,+,- at each trimer) for J' /J > 1. One can also apprcci- 
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describing a modulated local magnetization with period 
three. As shown in Fig. [11] (lower panel), this parameter 
allows the identification of three different phases. Notice 
that the second order transitions, characterized by long 
range fluctuations, are not expected to show up clearly 
in finite size systems. Thus the results around J'/J < 1 
are compatible with our previous semiclassical analysis. 

As a separate issue, we have also tested the confidence 
of the first order block-spin perturbation results. To this 
end, we computed the spectrum of the Hamiltonian ob- 
tained in Eq. fl"2"]) , adapted to a finite size system of 
24 spins with periodic boundary conditions, through the 
Jordan- Wigncr mapping to spinless fermions (care has 
to be taken in imposing periodic (antipcriodic) bound- 
ary conditions for odd (even) fermion filling to the re- 
sulting tight-binding Hamiltonian). The spectrum thus 
obtained can be compared with exact Lanczos diagonal- 
ization of the full spin system at different values J'/J. 
As it is known, the convergence of this method to highly 
excited states becomes progressively slow. To avoid this 
problem we restarted the Lanczos procedure with an ini- 
tial random state chosen orthogonal to each of the eigen- 
levels previously found. In Fig. [12] we plot in the up- 
per panel the first excitations for J'/J = 0.05, showing 
very good agreement both in values and degeneracy of 
the energy levels. In the lower panel the same is plot- 
ted for J'/J = 0.2, showing important deviations of the 
XY picture from the exact results. This deviation is ex- 
pected because of the level crossing found numerically at 
J'/J = 0.18 in the 24 spins system. 



FIG. 11: Upper panel (a): local magnetization of the ground 
state of 24 spins as a function of the coupling parameter. 
Squares (circles) correspond to middle 1 (ending 2, 3) sites in 
trimers. Panel (b) displays respectively these results for 18 
spins (here, the uniform phase is absent). Lower panel (c): 
order parameter in Eq. (I22|l describing modulated magneti- 
zations for 18 (circles) and 24 (triangles) spins. From left to 
right it exhibits uniform magnetization, quantum plateau and 
up-up-down ground states. 



ate a sudden magnetization change at the level crossing 
point J'/J = 0.18; the ground state magnetization pro- 
file corresponds to uniform magnetization (S z « 0.166 
at each site) for < J'/J < 0.18 and to the quantum 
plateau for 0.18 < J'/J. For the 18 spins system the 
magnetization pattern is similar, except for the smooth 
behavior between the uniform magnetization and quan- 
tum plateau phases; this corresponds to the above men- 
tioned possibility of the transition being BKT-like in- 
stead of first order. 

The numerical results can be summarized by an order 



V. SUMMARY AND CONCLUSIONS 

In the present note we have analyzed quantum phase 
transitions in zig-zag antiferromagnetic spin S = \ lad- 
ders at M = 1/3 driven by a trimerized modulation, as 
could be produced by adiabatic lattice deformations in a 
spin-Pcicrls like transition. 

Far from being stable, the triple degeneracy of the 
1/3 magnetization plateau in homogeneous ladders is 
lifted according to a triple sine-Gordon mechanism, giv- 
ing place to several magnetic phases separated by first 
and second order transitions. Numerical diagonalization 
results in finite systems support our conclusions. A fur- 
ther transition separates the weak modulation regime 
from the strong frustration regime. From weak trimer 
coupling to beyond the homogeneous zig-zag point, the 
phase diagram is schematically shown in Fig. [2] a phase 
with uniform magnetization described by non magnetic 
gapless chiral degrees of freedom, then a transition to 
a non degenerate quantum plateau state, then a second 
order transition, in the Ising universality class, to a two- 
fold degenerate state and finally a first order transition at 
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the homogeneous point to a non degenerate up-up-down 
ground state. The nature of the transition between the 
chiral gaplcss phase and the quantum plateau is not re- 
solved by our numerical data; we leave a blank in the 
phase diagram as an open question on this issue. 



Motivated by the spin-Peierls transition usually stud- 
ied at zero magnetization, it is interesting to investi- 
gate whether a lattice deformation like the one analyzed 
in this work could take place at low temperatures in 
magnetized systems because of the competition between 
magnetic energy gain and elastic energy loss. We hope 
that the knowledge of the ground state structure and 
non magnetic excitations presented here will be useful 
in such investigation. Related results will be published 
elsewhere. 
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FIG. 12: Comparison of non magnetic excitations \ n above 
the M = 1/3 plateau, obtained on a finite system from the 
first order block-spin perturbation and from exact diagonal- 
ization: (a) J'/J = 0.05, (b) J'/J = 0.2. The horizontal and 
vertical axes are the excitation level (counting degeneracy) 
and An — Ao respectively. 
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